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Abstract 

The Electron Localization Function (ELF) - as pro¬ 
posed originally by Becke and Edgecombe - has been 
widely adopted as a descriptor of atomic shells and cova¬ 
lent bonds. The ELF takes into account the antisymme¬ 
try of Fermions but it neglects the multi-reference charac¬ 
ter of a truly interacting many-electron state. Electron- 
electron interactions induce, schematically, different kind 
of correlations: non-dynamical correlations mostly affect 
stretched molecules and strongly correlated systems; dy¬ 
namical correlations dominate in weakly correlated sys¬ 
tems. Here, within an affordable computational effort, 
we estimate the effects of same-spin dynamical correla¬ 
tions on the electron localization by means of a simple 
modification of the ELF. 


I. INTRODUCTION 


At the core of quantum chemistry, there is the concept 
of chemical bond. Since early times, there have been 
ongoing efforts to set up effective descriptors of bonding 
in molecules and solids. Pair of electrons of opposite 
spins were proposed as the basis of bonding in the Lewis 
theory^. More recently the idea emerged that bonding 
and other structure, such as atomic shells, should reflect 
the “localization” of electrons in the system. 

An effective way to estimate the degree of localization 
is to exploit the fact that, if the position of one electron 
is fixed at a given point in space, it must be less likely 
to find a second electron around the same position. If 
the pair is made of electrons with parallel spin states, 
the localization should be related to the properties of 
the Fermi hole function which reflects the effect of Pauli 
exchange repulsion^. 

Becke and Edgecombe have followed these lines of 
thoughts and have analyzed a state of the form of a 
single Slater determinant^. In this case, the conditional 
distribution ''2) of having a cr-spin electron at 

position r2 given that there is a cr-spin electron at the 
reference position ri is given by 


-Pcond(ri, r2) = Pa{r2) + K{ri,r2) (1) 


(we employ a simplified notation with a single spin index 
as the pair of electrons is restricted to parallel spin state), 
where the exchange-hole (x-hole) function 


|Pi(i~Tr 2 )P 

Pa{ri) 


is expressed through the spin-dependent one-body re¬ 
duced density matrix 


Pi(ri,i'2) = Eir/>*^(ri)V'ia(r2) • (3) 

Here, the summation is restricted to occupied cr-state 
single-particle orbitals, ifiair), and the particle density 
is readily obtained as Pcr(r) = pf (r, r). 

Among the fundamental features of the x-hole func¬ 
tion, there are: 

• its normalization 

y d\2 /ix(ri,r2) =-1 , (4) 

i.e, it accounts for —1 electron; 

• its average extent 




\K{v,u)\ 


du 


(5) 


where 


/i^(r,M) := J dnhZir,r + u) (6) 

is the spherical average of the x-hole function de¬ 
termined around the reference position; 

• its short-range behavior for small u 

^x(j',u) = -Pair) - i [VVa(i-) - ^Dair)] 

+ ••■ (7) 


where the first term is also known as the on-top 
x-hole and the coefficient of the second term is the 
curvature of the x-hole. In detail 


Dair) 


Ta{r) 


1 (Vp.(r))' 

4 Pair) 


( 8 ) 


with 


Tair) := T;i\Vipiair)\'^ (9) 

being the (double) of the (positively defined) ki¬ 
netic energy density. 

Similarly to Eq. Q, 

Pair, u) = Pair) + ^V'^pair)u'^ H- (10) 
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and insertion of both Eq. © and Eq. OT in Eq. o 
yields 

Pcond{^,u) = ^Da{r)u^ + ■ ■ ■ . ( 11 ) 

Note, Eq. (ED is a key expression that we will later mod¬ 
ify [see Eq. (IT^ ] to take into account correlations beyond 
exchange effects to some extent. 

For an inhomogenous system, Du(r) exhibits a non¬ 
trivial dependence on the reference position r. It turns 
out useful to make a comparison with respect to a local 
uniform reference with same density as the actual system. 
Becke and Edgecombe introduced the quantity 

X.(v)=D,lv)/DT^{v) ( 12 ) 


Finally, we point out that the study of electron local¬ 
ization has been fundamental for the understanding of 
how to capture non-trivial non-localities in the exchange- 
correlation functional by means of simple approxima¬ 
tions. This and related topics continuously attract at¬ 
tention and novel insights have been discussed in recent 
contributions^^— . The analysis presented here adds use¬ 
ful information. 

The work is organized as follows: in Sec. m we present 
a modified ELF; in Sec. uni we analyze results of appli¬ 
cations to atoms, molecules, and Jellium; in the Sec. m 
we discuss conclusions and outlooks. 


II. MODIFICATION OF THE ELF 


where 


= (13) 

from which they defined the Electron Localization Func¬ 
tion (ELF) as 


■ («) 

The relevant information provided by the ELF is con¬ 
tained in its variations. The function varies from zero 
to one: for ELF = 0.5, the electrons are localized as 
in the uniform reference, otherwise they are more (less) 
ELF > 0.5 (ELF < 0.5) localized. The modulation of the 
ELF nicely reproduces the shell structure in atoms and 
well emphasizes covalent molecular bonds in molecules^. 

The ELF discussed above is built from a wavefunc- 
tion of non-interacting electrons. Works in the litera¬ 
ture have proposed generalizations of the ELF— at the 
full correlated level. Here, we find it useful to recon¬ 
sider this topic by distinguishing different type of correla¬ 
tions. While non-dynamical correlations affect stretched 
molecules and strongly correlated systems, dynamical 
correlations dominate in most of the molecules at typical 
equilibrium distances and weakly correlated systems. In 
the language of multi-reference methods, non-dynamical 
correlations are related to the mixing of configurations of 
low energy (which may be degenerate or almost degen¬ 
erate), while dynamical correlations have to do with the 
mixing of configurations at higher energy. Can one esti¬ 
mate the effects of dynamical correlations on the electron 
localization? Can these estimates be made within a well- 
affordable computational approach? As we illustrate in 
the next sections, these objectives can be achieved by ex¬ 
ploiting the ability of density functional approximations 
to capture semilocal correlations that are of dynamical 
type. The approach introduced here neglects opposite- 
spin contributions, not because they are expected to be 
unimportant but because we are concerned with an esti¬ 
mation implying a simple and direct modification of the 
original ELF. 


In an interacting system, we must account for the fact 
that electrons tend to avoid each other not only because 
of the antisymmetry requirement but also because of the 
electron-electron repulsion. As mentioned in the intro¬ 
duction, we shall focus on same-spin pairs, for which the 
conditional distribution is given by 


^cond(i'iT2) = Po{r2) + K{vi,T2) + K (fi. r2) (15) 


where (ri, r 2 ) describes many-body corrections beyond 
the effects included in the (Generalized) Kohn-Sham 
[(G)KS] state (notice that this definition differs from the 
coupling-constant averaged correlation hole). 

For a pair of electrons at small inter-particle distance, 
the interaction is strong yet the relative kinetic term is 
not suppressed. Thus, the pair in this configuration looks 
like a “hydrogen atom” but with a repulsive interaction 
and much lighter nucleousi^. Density functional theory 
(DFT) has successfully contributed to clarify that semilo¬ 
cal density functional approximations are very effective 
to capture such short-ranged correlations - that are, in 
fact, of dynamical type. 

Our choice is to resort to the model proposed (and 
explained in full details) by Beckei^ for the spherical av¬ 
erage of /i[f(ri,r 2 ) for r 2 ss ri. Becke’s final target was 
the DFT hole, whereas we are interested in the hole at 
full coupling strength: 


K{v,u) 


V? [u - Za{r)] Dajr) 

6 (1 -I- Zcr(r)/2) 


F(7„(r)u) , 


(16) 


where 


z,{r) := 2ci?5(r) , (17) 

has the meaning of a “correlation length” (it sets the 
shortest inter-particle distance at which the correlation 
hole vanishes at each reference position). In the follow¬ 
ing, we shall employ c = 0.88 as determined by Becke 
by fitting DFT correlation energies. The stability of the 
results with respect to relatively large variation of c will 
be demonstrated in the next section. F{x) is a damp¬ 
ing function with quadratic small-a: behavior {F(x) ^ 
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1 — ax^ + ...). 7 CT(r) will not enter in our final expression 
but it is instructive to appreciate that, it is determined 
by enforcing proper zero-normalization condition of the 
correlation hole and as a result, 7cr(r) ss l/za{r). For 
small u, Eq. (HU), reads 


hc(r,u) 


z„(r)D„(r) 2 

6 (1 -1- z^(r)/2)“ 


(18) 


The model in Eq. (HSl) has several appealing features: it 
sets a natural correlation length; it contains zero electron 
as the corresponding exchange-hole accounts for —1 elec¬ 
trons; it gives vanishing correlation energy for one elec¬ 
tron systems; and its short-range behavior is expected to 
be realistic particularly for inhomogenous systems. 

Upon insertion of Eq. (fT51) . Eq. (fTUl) . and Eq. ([7]) into 
(the spherical average around ri = r of) Eq. (fT^ . we get 


Tond(uw) « I <! D^{r) 


1 - 


Zair) 


2(l + z.(r)/2)J 


(19) 

from which we define a “modified” ELF (mELF) as fol¬ 
lows 

mELF(r) := --(20) 

1 + X?(r) 


where 


xUr) = 


X<T(r) 

I 

2 

X<t( 

r) 


Za{r) 


2(l + z^ir)/2) 


1 -b Zair)/2 ■ 


( 21 ) 


Notice that the numerator of Eq. (ED includes the ref¬ 
erence to the uniform gas as originally defined by Becke 
and Edgecombe. 

Finally, also notice that Eq. (17(11) - through Eq. (ITTl) . 
ED and Eq. m - entails the calculation of the Slater 
potential U^{r) 





( 22 ) 


i.e., the Slater potential is the potential due to the 
exchange-hole function. Features and implications of the 
modified ELF are explored in Sec. uni 


A. Further simplifications 

Since an evaluation of the Slater potential, U^{r), may 
not be readily available in some DFT codes, we remind 
that this quantity can be approximated in terms of the 
Becke-Russel (BR) model Ugj^gg ■ The BR poten¬ 

tial is appealing for finite systems because it entails a 
model for the x-hole which is exact for finite 1-particle 
system and reproduces exactly the shot-range behavior of 
/ix(r, u) in Eq. ([T]) of any fV-electron system. One should 


keep in mind that, by construction, it produces a poten¬ 
tial of a “localized” hole function. Within our scope, the 
BR model can be straightforwardly applied using orbitals 
obtained from a converged (G)KS calculation. 

In the case of molecules, where the exchange-hole can 
delocalize over multi centers, the option to adopt the full 
Slater potential is available from direct post processing 
of the orbitals obtained from a converged (G)KS calcu¬ 
lation. 

For applications in extended metallic-like systems, we 
may consider the LDA expression for the Slater potential 


xrLDA,(T 


-3 




(23) 


This expression is exact for Jellium. 


B. Current-carrying and time-dependent states 

In the expressions above, we have tacitly assumed real¬ 
valued orbitals. For current-carrying states in the case of 
ground state degeneracy—^— and time-dependent pro- 
cessesi^, complex-valued orbitals must be considered in¬ 
stead. The proper expressions are obtained with the 

substitution r(r) T(r) — > where jp(r) is 

the paramagnetic current. Of course, a time-dependence 
(r) —>• (r, t) must also be admitted for the time-dependent 
case. In the reminder of this work, we shall restrict our¬ 
selves to systems with closed-shell ground states. 


III. APPLICATIONS AND FURTHER 
ANALYSES 

In this section, we apply Eq. m and Eq. (EOD to sev¬ 
eral systems in order to quantify their differences. 


A. Atoms 

As a first system, we consider the Argon atom. Fig. [1] 
reports the ELF and its modified version (i.e, mELF). 
Qualitatively, the resulting pictures of the shell structure 
exhibits strong similarities. The included short-ranged 
correlations do not change the position of the maxima 
and minima. However, we can appreciate quantitative 
differences: minima are less shallow and maxima higher. 
Localization is enhanced and it increases and extents as 
moving outwardly from the center of the atom. 

Fig. [T] was obtained by employing the BR model for 
the Slater potential. In order to make the illustration 
self-contained, we report the BR and Slater potential for 
the Ar atom in Fig. [2] 

Finally, the dependence of the results on the values of 
the parameter c of Eq. ([T7|) is studied in Fig. [H It is ap¬ 
parent that variations of about ±20% do not change the 
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FIG. 1: (X) Electron Localization Function in its original 
form, Eq. m, and (X+C) in its modified version, Eq. (I20II . 
computed for the Argon atom. Kohn-Sham orbitals ob¬ 
tained with the APE code^, using a LDA for the exchange- 
correlation energy functional, are used as input to the two ex¬ 
pressions. The Slater potential, which enters Eq. ([20]) through 
the correlation length Za- [see also Eq. (|17l) and Eq. (1221) ]. is 
evaluated within the BR approximation (see also Fig. [2J. 



r (a.u.) 


FIG. 2: Dots show the Slater potential for the Argon atom 
evaluated within the OEP-KLI method^ as implemented in 
the OGTOPUS code^ (a grid with equally spaced points is 
employed). The solid line shows the same quantity evaluated 
within the BR approximation (this was then employed in the 
calculation of the modified ELF reported in Fig. [T|. The BR 
potential is evaluated on quantities produced with the APE 
code as explained in the caption of Fig. [T] It is apparent 
that, in the region of interest, the two potentials do not differ 
significantly. 

results significantly. Hence, we do not find any evidence 
of the necessity to optimize the value of c further. 


B. Molecules 

In order to assess the effect on the ELF of dynamic cor¬ 
relation, as introduced in this work, we have compared 
the usual ELF with the mELF for a few small hydrocar¬ 
bons with different CC bond orders (Ethane, Ethene and 


FIG. 3: Demonstration of the stability of the results reported 
in Fig.[T]for the modified ELF [ Eq. (1201) ]. Here, the parameter 
c - entering the final expression through Eq. 113 and Eq. 123 
- is varied ±20% around its prescribed valued. 

Ethyne). 

We have verified that the pictures obtained with the 
standard ELF is essentially unchanged with respect to 
what reported, for instance, in Ref.—. The corrections in¬ 
troduced by mELF are only minor in these cases. More 
noticeable differences may be found in molecules involv¬ 
ing heavier atoms. 

Let us consider the case of Iodine molecule in Fig. S) 
The qualitative behavior of mELF follows that of the 
ELF with larger values in the bonding region. We have 
verified that, when the two atoms are brought closer and 
closer, the relative minimum at the midpoint of both ELF 
and mELF is turned into a maximum when the bonding 
distance is reached (approximately 2.7 A). At this point, 
the electrons pairing from the two 5p atomic orbitals be¬ 
comes a TT type molecular orbital. For all the molecules 
considered in this section, we have input the mELF with 
the Slater potential evaluated with converged LDA KS- 
DFT results. 


C. Jellium 

The application of Eq. (uni) to Jellium (an extended 
uniform gas of interacting electrons, whose charge is bal¬ 
anced by a smeared positive background) gives us the op¬ 
portunity to put forward some explorative speculations. 

For this system, the original ELF is a constant inde¬ 
pendent of the particle density (the values of the constant 
being 0.5). Instead, our mELF gives a constant that de¬ 
pends on the values of the particle density. Eq. (II3 to¬ 
gether with Eq. (l22l) and Eq. (l23l) show that the correla¬ 
tion length is in inverse relation with the particle density 
z^ir) - pa^^^(r). 

For Jellium, the low-density limit corresponds to the 
strongly interacting regime. If we require the KS wave 
functions to be plane waves, we find that the correlation 
length diverges Za-(r) —>■ ±00 and thus mELF ^ I. Since 
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FIG. 4: (X) Electron Localization Function in its original 
form, Eq. m, and (X+C) in its modified version, Eq. (I20II . 
for the Iodine molecule. The calculation was carried out with 
the OCTOPUS code using a LDA for the exchange-correlation 
energy functional. In (X-l-C), we have employed the Slater 
potential obtained post-LDA. We employed Troullier-Martins 
norm-conserving pseudo potentials (as the focus here is on 
the state of the valence electrons). 

it is believed that, in the same limit, the electrons in 
Jellium tend to localize into a crystalline structure (i.e., 
the Wigner’s crystal)^, we may say that, mELF captures 
this instability. 

In the opposite limit, electrons in Jellium get weakly 


interacting. Correspondingly Za{r) —>■ 0 and thus 
mELF —)• 0.5. We may say that, electrons fully de¬ 
localize in order to get packed together to develop the 
high-density limit. 


IV. CONCLUSIONS AND OUTLOOKS 

We have shown how to account for the effects of 
same-spin dynamical correlations on the electron lo¬ 
calization within a model which requires an affordable 
computational effort. We could visually asses that the 
effects of these correlations are somewhat unimportant 
in small organic molecules. They imply some noticeable 
effects for atomic shells and bonds involving relatively 
heavier atoms. However, qualitative features such as 
the positions of the shells and the bond character 
obtained at the exchange-only level are unchanged. 
As a novel interesting feature, the proposed modified 
electron localization function appears to connect the 
degree of localization in Jellium to different interaction 
regimes. For the future, it is appealing to attempt to 
include information on the opposite-spin channels and 
of non-dynamical correlation effects. 
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